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Abstract 

The techniques of integration by parts and differential reduction differ in the 
counting of master integrals. This is illustrated using as an example the two- loop 
sunset diagram with on-shell kinematics. A new algebraic relation between the mas- 
ter integrals of the two-loop sunset diagram that does not follow from the standard 
integration- by-parts technique is found. 
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Within dimensional regularization [T], the integration-by-parts (IBP) technique 
[2] (for a recent review, see Ref. [3]) is one of the most powerful tools for evaluating 
multiloop Feynman diagrams. However, the question how to construct algebraical, 
geometrical, or topological criteria for the (ir)reducibility of Feynman diagrams in 
general is still an open one At present, only the explicit solution of recurrence 
relations (see, for example, Ref. [5]), or the application of the Laporta algorithm 
[6], as implemented in computer programs [7J, provide an answer on this question. 

In a series of recent publications [8 9 l O llj . we addressed the problem of count- 
ing nontrivial master integrals with the help of hypergeometric representations for 
Feynman diagrams and the differential-reduction algorithm [12]. By detailed anal- 
ysis, we thus derived the following empirical criteria for the hypergeometric repre- 
sentation of Feynman diagrams: If a Feynman diagram is expressible as a linear 
combination of Horn- type hypergeometric functions with rational coefficients, then 
(i) each hypergeometric function has the same number of basis elements in the 
framework of differential reduction [l2] ; and (ii) the number of nontrivial master in- 
tegrals is equal to the number of basis elements of the hypergeometric functions (up 
to the modules of Feynman integrals expressible in terms of products of algebraic 
functions and F functions). 

These criteria were conjectured on the basis of the analysis of a variety of par- 
ticular examples [8|9|10j . A rigorous mathematical proof is still lacking. The aim of 
this Letter is to present a first example where the application of our criteria allow 
us to find an additional algebraic relation between master integrals which does not 
follow from the standard IBP technique. In fact, Tarasov's algorithm [5] for the re- 
duction of two- loop propagator diagrams, as implemented in the computer packages 
of Ref. [13] , and Laporta's algorithm [6] , as implemented in the computer package 
of Ref. [7], fail to produce this relation. 

Let us consider the two-loop self-energy sunset-type diagram J012 with on-shell 
kinematics, defined as 



where n = 4 — 2e is the dimensionality of space time (see Fig. [T]) . Such a diagram 
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Figure 1: Two-loop self-energy sunset-type diagram Joi2- 
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contributes to the pole mass of the top quark [T3]. The hypergeometric representa- 
tion of this diagrams was presented in Eq. (3.12) of Ref. [H]. We reproduce it here 
for completeness: 
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In the framework of differential reduction [9], the first hypergeometric function in 
Eq. ([2]) is expressible in terms of a Gauss hypergeometric function and a rational 
function Ri{z) of z = 4m^/M^, whereas the second one is expressible in terms 
of a 3F2 function with one unit upper parameter and a rational function R2{z). 
Schematically, we have 
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where = zd/dz, the short-hand notation (1, 9) stands for {Pi{z) + P2{z)9) , with Pi 
being rational functions (see Eqs. (17) and (20) in Ref. [9j), and Ij (i = 1, . . . , 4) are 
integers. According to our criteria, there are two master integrals for this diagram 
that are not expressible in terms of F functions. However, as shown by Tarasov in 
Ref. [5], solving the standard IBP relations [2] yields three master integrals of this 
type, namely Joi2(l,l,l), '^012(1, 2, 1), and Joi2(l,l,2), which is confirmed with 
the help of the computer packages of Refs. [7|13j . Consequently, either our criteria 
are wrong or there exists an algebraic relation between the integrals Joi2(l,l,l)j 
Joi2(l, 2, 1), and ^012(1, 1, 2), possibly including some algebraic functions depending 
on z and products of V functions. 

To find this relation, let us explore Eq. ([2]) and present the master integrals in 
the following form: 



Xiii = Axi + Byi , X121 = Ax2 + By2 , X112 = Ax^ + By^ , 



where 
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xi = 2F1 {^,a;b;z) , X2 
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with 



a = 3 — n , "~2^' "^"2 — — , a = n — 1 , (6) 

and pFp-i{a; b; z) being a hypergeometric function. Notice that the 4-F3 functions in 
Eq. [2] collapse to 3-F2 and 2F1 functions according to Criterion I defined in Section 2.4 
of Ref. [9j because upper indices exceed lower ones by integers. Using the relations 

apFq{a + l,A;B-z) = {6 + a) pFg{a, A; B; z) , 
(6 - 1) pF,{A; b-l,B;z) = {e + b-l) pFg{A; b, B; z) , (7) 

it is easy to obtain 

{3n — 8)xi + ZX2 + 2x3 = n — 2 , 

(3n - 8)yi + zy2 + 2ys = 0. (8) 

This is equivalent to the following relation between master integrals: 

(3n-8) Joi2(l, 1, 1) + 4m2 Joi2(l, 2, 1) + 2M2 Joi2(l, 1, 2) 

= 2(M2r-3rg-l)r(3-n)r(2-^) . (9) 

Eq. ([9]) does not only represent a relation between Feynman diagrams, but also 
a relation between hypergeometric functions pFg, which is newly derived from the 
differential reduction technique. 
For m? = M^, we have 



(3n - 8)Joii(l, 1, 1) + 6mVoii(l, 1, 2) = 2(771^)"-^? (3-n) P [2-- 

(10) 

The integrals Jqh (1,1,1) and Joii(l,l,2) with m? = are master integrals of 
the package 0N-SHELL2 [15j. Eq. ^ coincides with Eq. (4.45) of Ref. [I6], where it 
was derived by studying the analytical coefficients of the higher-order e expansion 
of diagram Joii(l,2,2) performed in Ref. [17j . 

According to Ref. [18], the last terms in Eqs. ([9]) and (llOp may be identified 
with a two-loop bubble diagram with one massive line divided by M^(n — 2), so 
that Eq. ([9]) is an algebraic relation between four two-loop diagrams. However, that 
two-loop bubble diagram, having two massless lines, cannot be obtained from the 
original two-loop sunset diagram with two massive lines by the contraction of lines. 
Consequently, Eq. ([9|) cannot be derived from standard IBP relations for two-loop 
sunset diagrams. 

At this point, a comment on the structural difference between the differential- 
reduction and IBP approaches seems appropriate. Eq. ^ represents a linear relation 
between master integrals that contains an inhomogenious term involving P functions 
on the right-hand side. However, it is clear from the structure of the IBP relations 
that the IBP method can only result in relations between integrals in which the 
prefactors are rational functions, while P functions do not appear in IBP relations. 
This is another reason why Eq. ^ cannot be derived using the IBP method. 
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In conclusion, counting the number of irreducible master integrals for the dia- 
gram defined by Eq. ([T]) using our criteria on the one hand and the explicit analytical 
solution |5|6j of the standard integration- by-part relations [2] , as implemented in the 
widely used program packages [7,13j, on the other hand, we encountered a mismatch. 
This motivated us to established a new algebraic relation, namely Eq. ([9]). This is 
a first example where our criteria allow us to expose a new relation between master 
integrals of the standard integration-by-part technique. 

Finally, we wish to mention that the final result of Ref. [Tlj does not depend on 
the number of master integrals. In fact, all master integrals are expressible in terms 
of hyper geometric functions, and the analytical coefficients of the e expansions of 
the latter were constructed in Ref. |19j . 
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